Beating the standard quantum limit for force sensing 
with a coupled two-mode optomechanical system 
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Optomechanics allows the transduction of weak forces to optical fields, with many efforts ap- 
proaching the standard quantum limit. We consider force-sensing using a mirror-in-the- middle setup 
and use two coupled cavity modes originated from normal mode splitting for separating pump and 
probe fields. We find that this two-mode model can be reduced to an effective single-mode model, 
if we drive the pump mode strongly and detect the signal from the weak probe mode. The optimal 
force detection sensitivity at zero frequency (DC) is calculated and we show that one can beat the 
standard quantum limit by driving the cavity close to instability. The best sensitivity achievable is 
limited by mechanical thermal noise and by optical losses. We also find that the bandwidth where 
optimal sensitivity is maintained is proportional to the cavity damping in the resolved sideband 
regime. Finally, the squeezing spectrum of the output signal is calculated, and it shows almost 
perfect squeezing at DC is possible by using a high quality factor and low thermal phonon-number 
mechanical oscillator. 



Dramatic progress in coupling mechanics to light [1-4] 
suggests that such devices may be used in a wide vari- 
ety of settings to explore quantum effects in macroscopic 
systems. Furthermore, such systems can be exquisitely 
sensitive to small perturbations, such as forces induced 
either by acceleration as in accelerometer [5] or by, e.g., 
coupling to surfaces or fields as in atomic force mi- 
croscopy [6]. For such force measurements, a high quality 
factor (Q) mechanical oscillator acts as a test mass, trans- 
ducing a force into a time-dependent displacement of the 
oscillator [7, 8]. By using intcrfcromctric techniques to 
monitor the position of the oscillator, one can infer the 
force via optical signals. However, the radiation pres- 
sure coupling between the mechanical mode and optical 
mode has three consequences: photon shot noise, quan- 
tum backaction and dynamical backaction [1, 9]. The 
dynamical backaction modifies the oscillator dynamics 
[10] and makes laser cooling [11, 12] or amplification of 
phonons [13] in the mechanical system possible. Pho- 
ton shot noise and quantum backaction, the former de- 
creases with increasing input laser power while the latter 
increases with increasing input laser power, introduce two 
sources of noise on the displacement readout of the oscil- 
lator motion. An optimal compromise between these two 
noise sources leads to the standard quantum limit (SQL) 
in force sensing [7]. 

The SQL, however, is itself not a fundamental limit. 
By using squeezed states of light [14], employing quan- 
tum nondcmolition (QND) measurement [15], or by cav- 
ity detuning [16], the SQL can be surpassed. Here we 
show that in a coupled two-mode optomechanical system, 
if we drive it appropriately, the interaction between cav- 
ity photons and the mechanical oscillator will generate 
squeezed states of the output light. Measuring an appro- 
priate quadrature of the output light field, we would get 
fewer fluctuations than that of the vacuum state, which 



makes it possible to detect weak forces far below the 
SQL. Furthermore, since we pump and probe different 
resonant optical modes, the effective optomechanical cou- 
pling is enhanced, and thus the pump power requirement 
for achieving the best sensitivity is lowered substantially. 

We consider a high finesse Fabry-Perot cavity with a 
dielectric mirror in the middle [17, 18] [Fig. 1(a)]. The 
end mirrors of the FP cavity are fixed, while the middle 
mirror can vibrate along the optical axis of the cavity at 
a mechanical frequency u> m with effective mass m. With 
the presence of the middle mirror, the FP cavity is di- 
vided into two sub-cavities, denoted by left (L) cavity 
and right (R) cavities. The middle mirror has a nonzero 
transmission, which allows the exchange of light between 
these two subcavities and thus leads to an effective cou- 
pling between the left and right cavity modes [19]. Fur- 
ther, the coupling will shift the resonant frequencies of 
two coupled cavity modes and leads to the so-called nor- 
mal mode splitting effect [20, 21]. 

Following [18, 22], the normal mode splitting in the 
presence of the middle mirror can be calculated by assum- 
ing a transfer matrix with a high reflectivity ra ■ For sim- 
plicity, we assuming the middle membrane to be exactly 
at the middle point of the FP cavity initially, dividing the 
cavity into two subcavities with the same length L, and 
the normal mode splitting when the middle mirror is at 

c 

a new position x is given by il = — arccos(|rd| cos(2fccc)) 

L 

[Fig. 1(b)], where is complex amplitude reflectivity of 
the middle mirror, k is the wave vector of the incoming 
field and x is the displacement from the middle point of 
the FP cavity. 

The hamiltonian of the cavity fields in this three- mirror 
system is 
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FIG. 1. (Color online), (a) A symmetric mirror-in-the-middle 
optomechanical system comprising a Fabry-Perot cavity of 
length 2L with a high reflectivity mirror mounted in the mid- 
dle and coupled to a mechanical oscillator. Displacements of 
the middle mirror (via mechanical oscillations) couple the two 
normal modes a (red) and b (blue) as the left-right symmetry 
is broken, (b) Normal mode frequencies ui cav ,± (green) as a 
function of middle mirror displacement x. (c) Transmission 
spectrum of the three-mirror cavity, showing pairs of normal 
modes. We drive mode a strongly (long red line) and detect 
mode b (short blue line). 



where oj c is the resonance frequency of the subcavities 
with the middle mirror exactly in the middle (x = 0), / = 
{d(jJ C av/dx)\ x=Q is the shift of cavity resonant frequency 
per unit length evaluated at x = 0, and at (or) is the 
left (right) cavity mode annihilation operator. Here only 
the linear order frequency shift is considered, since the 
displacement x is much smaller than the cavity length L. 
The last term describes the coupling between left/right 
cavity modes with strength g. In the high reflectivity 
limit (|r<j| — > 1), g — \td\c/2L and / = —U) C JL, where 
td is the amplitude transmission coefficient and c is the 
speed of light in vacuum. 

We look only in a narrow spectral range around a 
nominal pair of normal modes a = [cll + afl)/v2 and 
b = (ol — aR)/y/2. We also drive mode a strongly at 
frequency col and move to the rotating frame with re- 
spect to the pump laser [23]. The cavity field is cou- 
pled to fields outside the cavity through the ends mir- 
rors, while we assume the mechanical oscillator is cou- 
pled to a thermal bath through clamping losses. The 
classical and quantum fluctuations of the environmental 
degrees of freedom will introduce damping to the cavity 
field and mechanical oscillator [24, 25], as required by 



the fluctuation-dissipation theorem [26]. In the Marko- 
vian approximation, the Hciscnbcrg-Langevin equations 
for mechanical and optical degrees of freedom are, in the 
high temperature limit, as follows: 
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In the equations above, 7 is the damping of the me- 
chanical oscillator, k is the damping of the cavity, and 
A c = lo c — lul is the cavity detuning. E is the pump 
strength which is related to input laser power Pj n and 
cavity damping k by \E\ = ^/P m n/huj L . a ln , b in are 
the vacuum fluctuations of the two cavity modes. Fi n is 
the force acting on the oscillator, and it has two parts: 
an external force F ext acting on the oscillator, which is 
also the force to be detected [Fig. 1(a)]; the Brownian 
stochastic force, or thermal fluctuating force F t h, which 
leads to damping of the oscillator. 

We then find the steady state of the system to be (x) = 
0, (p) = 0, (6) = p = 0, (o) = a = E/[i(A c - g) + «], 
which is the only stable solution at low power. Following 
[27, 28], we study the stability of this solution by apply- 
ing the Routh-Hurwitz criterion [29], which at positive 
effective detuning A = A c + g > is simplified to a 
constraint on pump strength, 



< »i^L 2 (k 2 + A 2 )/2hAuj 2 c = a 
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Writing x, p, a, 6 in terms of steady state value and fluctu- 
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ation, and neglecting terms of order — bx, — ax, ab, 

L L L 
then the fluctuation of mode a decouple, and the equa- 
tions of motion reduce to, 

x = p/m, (4a) 
p= -rnuj 2 m x~"/p + hG(b + b r )+F in , (4b) 
6 = -iAb - Kb + iGx + V2~Kb in , (4c) 

where we choose a to be real and define G = ou c a/L. 
For convenience, we remove the tilde of the fluctuating 
variables. We find that the pumped, coupled two-mode 
model is reduced to an effective single-mode model [30], 
where the cavity consist of a fixed partial transmitting 
mirror and movable perfect reflecting mirror and the op- 
tomechanical coupling strength is determined by a, the 
strength of pump. 

We then define the quadratures of mode b as X = 
(6+6t)/\/2, Y = (b — tf)/iy/2, move to frequency domain 
by Fourier transform, and solve a set of linear equations. 
We find 



x (u) = x(u)F in (oj) + x(w) 
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where x(w) is the susceptibility of the optomechanical 
system to force, 
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From the expression above, we immediately identify an 
effective, frequency-dependent mechanical resonant fre- 
quency uj' m and an effective damping 7' which are shifted 
from the original ones. The shift in resonant frequency is 
the so called "optical spring" effect [10], while the shift 
in damping leads to cooling or heating of the oscillator, 
depending on the sign of detuning [1, 2]. The cavity field 
fluctuations enter the equation of motion for oscillator 
Eq. (5) as an additional fluctuations force, which is iden- 
tified as the shot noise fluctuations of radiation pressure 
force. 

Within the input-output formalism [25], the output 
field quadratures we measure are related to the field 
quadratures inside the cavity by X out (uj) = ^/2kX(uj) — 
X in (uj) and Y out (uj) = ^/2kY(w) — Y in (uj). We consider 
a homodyne measurement of the signal [25, 31] 



S(lu) = sin 9 X out (lu) + cos QY 0Ut {u) 
= XF(w)F in (u)) + xx{u)X in (uj) 
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Xy{u)Y in (uj) 



where is an experimentally adjustable phase, which de- 
termines the measured quadrature. Here the signal is 
written in terms of three inputs (force input Fi n , ampli- 
tude fluctuations Xi n , phase fluctuations Yi„) and corre- 
sponding susceptibilities. The force and field susceptibil- 
ities are: 
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To calculate the sensitivity to the external force F ext , 
we define the following quantity 



F(u) = 



dS(oj)/dF e: 



(9) 



and the square of force detection sensitivity is given 
by the power spectral density of F(uj): t){uj) = 
J cLu' (F(lo)F(lo')) . The vacuum radiation input noise bi n 



is delta correlated and the thermal fluctuating force is ap- 
proximated as white noise thus is also delta correlated, 
so we have 



rj(uj) = 2m~fk B T + - 



Xf(w) 



(10) 



The first term is the thermal white noise. There is a 
term proportional to G 2 which comes from the random 
back action force, and also a term proportional to 1/G 2 
which is the phase noise related to position measurement 
imprecision. 

We now focus on the DC (u> = 0) force sensing regime. 
To get the best sensitivity, we first optimize the func- 
tion rj(u) = 0) for a and then optimize for 9. We find 
that for the optimal pump strength a 2 , the second term 
goes to zero as Asin# + kcos9 — > 0, which corresponds 
to a backaction free point. At first sight, it seems that 
Xf approaches zero as Asin6> + kcos9 — > 0, and the 
sensitivity diverges. However, if we choose the pump 
strength appropriately, then not only can the divergence 
at 6*o = — arctan(«;/A) be avoided, but also the sensi- 
tivity can achieve its optimal value. The optimal pump 
strength is given by 



oil 



1 



2k A sin 9 + k cos 9 



A V« 2 + A 2 
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where ao is the threshold pump strength defined in 
Eq. (3). At this point, the effective mechanical frequency 
uj' m — » as 9 — > 9 , so xf is still finite. We then try to 
find out the behavior of the sensitivity near the criti- 
cal angle #o- To ensure that the pump strength does 
not exceed the threshold value and that the effective 
mechanical frequency is positive, we let the angle 9 ap- 
proaches 9q from the positive side, that is 9 = 9q + 59, 
with < 69 <C 1. We note that, if 9 approaches 9q 
from the negative side, we can replace the minus sign in 
Eq. (11) with a plus sign, and the result will be similar. 
In the vicinity of 9q, the pump is approximated as 



(12) 



At the optimal pump strength and optimal angle, the 
total sensitivity is found: 



= 0) w 2m7&BT + hmuj 2 



A K 

4k + A 



(13) 



with the dimcnsionless parameter £ defined by £ = 
59 ■ 2k/ A. This result implies that thermal noise lim- 
ited detection can be achieved by choosing the critical 
angle 9q appropriately. The expense for achieving the 
best sensitivity is that we have to pump the system at 
a power close to the threshold value, thus increasing the 
possibility of destabilizing the system. The ultimate sen- 
sitivity for force detection is limited by how strong the 
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thermal noise is and how close we can pump the sys- 
tem near its instability point. The result also suggests 
that we could further improve the sensitivity by choosing 
A = 2k. Then we look at the input laser power, and find 
that if A c = g, the power is minimized. Along with the 
condition A = 2k, we have A c = g = k, and the optimal 
pump power is 



(14) 



where Q c = lo c / k is the quality factor of the cav- 
ity. Considering an optomechanical system with m = 
5.36 x 1CP 10 Kg, uj m = 2n x 130 kHz, pumping laser 
wavelength A = 1.55 /im, and cavity finesse of F = 20000 
[32], we find that for £ = 0, P op t = 0.816 mW, and the 
circulating power P cn - = P op tF/27r = 2.56 W. 

In practice, we need to understand the behavior of r)(w) 
at low but nonzero frequencies to determine the band- 
width for force detection. A full analysis is only possible 
numerically, so here we present a simple but illuminating 
approximation method. The argument is the following: 
at nonzero frequencies, in order to keep the 58 depen- 
dence in Eq. (13), we require \u> cos9\ < A s'md + k cos 9, 
which is equivalent to 



I* 

K 



(15) 



Thus the bandwidth is approximately [l + (A/k) 2 ] £k. 
Neglecting the thermal noise term, we plot the DC sen- 
sitivity for different values of £ in Fig. 2(a), and the sen- 
sitivity at low signal frequencies for different values of k, 
f, A/k in Fig. 2(b), 2(c), 2(d) respectively. We find that 
the bandwidth shown in the numerical result is in good 
agreement with Eq. (15). Thus we could use this formula 
to estimate the bandwidth at low frequencies for a given 
set of parameters. 

Finally, we study the squeezing spectrum of the output 
signal S(oj), which is given by S(uj) = / du' (S(ui), S(ui')) 
[25, 33]. Using the noise correlation relations, we have 

S(u) = 2m~fk B T\xF(u)\ 2 + - |%xM - ixv{u)\ 2 • (16) 

At DC, we can minimize it by choosing the optimal angle 
and pump strength, obtaining 



(17) 



where n t h = &bT '/hw m is the number of thermal phonon 
and Q m = uj m /"f is the quality factor of the mechani- 
cal oscillator. For A = 2k, S = 2(1 - C)n th /Q m + £ 2 . 
Nano-mechanical oscillators of high quality factor and 
low phonon number have now been fabricated by many 
groups, which makes it possible to reduce the thermal 




FIG. 2. (Color online), (a) Sensitivity at DC as a function 
of A/k for different values of £ (numerics shown in plot). 
At A/2/c, the minimum value is achieved. In the following 
figures, we choose u m — 2ty x 10 6 Hz, k = 0.2oj m , A = 2k and 
£ = 0.01 as the base values to plot the sensitivity as a function 
of frequency. The bandwidth is the frequency range where the 
best sensitivity is maintained, (b) Bandwidth dependence 
on k, where uj m is fixed and k varies and the ratio n/cj m is 
shown in plot. In the resolved sideband regime the bandwidth 
increases with k linearly, while in the unresolved regime, it 
does not increase much, (c) Bandwidth dependence on £, 
which is approximately linear, (d) Bandwidth dependence 
on A/k , where k is fixed and A varies. The dependence 
is approximately in the form yT~+ (A/k) 2 and suggests a 
narrower bandwidth than our estimate Eq. (15). 



noise term to a very small value. By driving the sys- 
tem near the threshold (£ — > 0), squeezing (5 < 1) in 
the signal we measure could be realized. This is the op- 
tomechanical analog of the squeezing in the output field 
from an optical parametric oscillator (OPO) [23]. Simi- 
larly, choosing 8 — > ir/2 + 8, we calculate the squeezing 
spectrum of a signal ir/2 out of phase and find that the 
optical noise term is proportional to l/£ 2 , which is con- 
sistent with the Hciscnberg uncertainty relation. 

In all the analysis above, however, we have not yet 
taken into account the quantum noise X' in (ui) introduced 
when measuring the field quadrature S(u>). Considering 
a measurement efficiency P < 1, then the actual signal 
we measure is S'(u) = \/PS(lu) + ^/T^ TX' m {uj). We 
find the modified sensitivity to be ?/(w) = t)(lj) + (1 — 
P)/2P\\f(uj)\ 2 . Assuming that the optimization at a 2 
and 8q is still valid, then the sensitivity at DC is rj'(0) = 
77(0) + ftmw 2 ra (l — P)/(l — £) • k/A, and the corresponding 
squeezing spectrum is 5" = S + (1 — -P)/2, in the limit 
P — > 1. Thus, a homodync measurement efficiency of 
99% will limit the squeezing to 23 dB (decibel). 

Our approach to squeezing for improved force sensing 
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may also have direct application in related topics, includ- 
ing atomic force microscopy, magnetic resonance force 
microscopy, and even in quantum transduction via me- 
chanics as recently suggested [34]. Furthermore, more 
complicated cavity mode structures, such as those of 
higher orbital angular momentum, may provide addi- 
tional methods for achieving this outcome in single-side 
cavities, as the fundamentals of our approach are simply 
having two well isolated, near-by cavity modes that both 
interact with the mechanical degree of freedom. 
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SUPPLEMENTAL INFORMATION 

NORMAL MODE SPLITTING IN THE HIGH REFLECTIVITY LIMIT 

The hamiltonian of the cavity fields in this three-mirror system is 

H opt = h(oj c + fx)a} L a L + h(ui c - fx)a R a R - Kg(a' L a R + a R a L ), (18) 
Diagonalizing this coupled left/right mode hamiltonian, we find two new eigcnfrcqucncies: 

u cav>± {x) = uj c ± y/px 2 + g 2 , (19) 

which correspond to the normal mode frequencies of the three-mirror system, and the difference between these two 
eigenfrequencies is exactly the normal mode splitting in Eq. (19), 

2 v / f 2 x 2 + g 2 = y arccos(|r d | cos(2fcx)). (20) 

Matching the left and right hand sides at a small displacement range, we find g = arccos(|rd|)c/2L and / = 
-y/\r d \ a,rcsm(\t d \) /\t d \uj c / L. In the high reflectivity limit (\r d \ -> 1, \t d \ -> 0), g = \t d \c/2L and / = -u c /L, 
which implies that the coupling constant is proportional to the transmission of the middle mirror and that the dis- 
persion of the cavity resonant frequency is linear. We have to mention that the full expression for the resonant 
frequency shift / is valid at both high and moderate reflectivity, as soon as the displacement x is small compared with 
the wavelength of the cavity field. But problem rises at moderate reflectivity, because the normal mode splitting 
increases as reflectivity decreases, and it approaches the free spectrum range of the large cavity with length 2L. To 
avoid the coupling between multiple cavity modes, we only consider the high reflectivity case, which has been realized 
experimentally [35]. 

To simplify the hamiltonian, we define two new modes a = (a_L + or)/\/2 and b = (ai — a R )/V2. If we drive mode 
a strongly at frequency ojl, we can move to the rotating frame with respect to the pump laser, then we have the 
hamiltonian of the whole system (vibrating mirror, cavity modes and pump): 

T) 2 X 

H = — h -rmo 2 x 2 + h(A c - g)<va + h(A c + gWb - tno c — (a% + Va) + ihE(a) - a), (21) 

2m 2 L 

where A c = uj c — ojl is the cavity detuning, and E is the pump strength which is related to input laser power Pi n and 
cavity damping n by \E\ = y Pin^/fv^L- 



SOLVE THE EQUATIONS OF MOTION 



To solve the reduced Heisenberg-Langevin equations, we define the quadratures of mode b as X = (b + f>t)/v2, 
Y = (b — b^)/iy/2, and move to frequency domain by Fourier transform, obtaining the following equations of motion: 



- iu!x(uj) = p(ui)/m, 

r % 2 / \ / \ V2fiaj r a . . . 
-iujp(uj) = -muj m x{Lj) - 7p(w) H X(u) + b in {ijj), 

-iojX(lo) = -kX(lj) + AF(w) + V2~KX in (uj), 
-iLuY(uj) = -kY(u) - AI(u) + ^ cOL x{uj) + V2~KY in (u}). 



L 



From this set of coupled linear equations, we find, 



¥(«,) 



1 



(k - iuj) 2 + A 2 
1 

(k - itu) 2 + A 2 



V2uj c aA 
L 



x(uj) + V2~k(k — iui)Xin(u)) + V2~KAYi n (uj) 



(k — iu))x(u>) + v 2k(k — iLo)Yi n (uj) — \J2tiAXi n (uj) 



(22a) 
(22b) 
(22c) 
(22d) 

(23a) 
(23b) 



At zero detuning A = 0, X(uj) is unchanged, while Y(io) is modulated by oscillator displacement x(uj) [1]. Measuring 
Y(lo) will give us the information about oscillator displacement. At finite detuning, however, both quadratures are 



7 



related to oscillator displacement. Putting the field quadratures back into the equation of motion for the oscillator, 
we have 



■[(« - iuj)X in (uj) + AY in (u)] 



(k — iu>) 2 + A 2 

with G = ui c a/ L and x( w ) is the susceptibility of the optomechanical system to force, 

2hG 2 A/m 



(k — ilo) 2 + A 5 



(24) 



(25) 



FORCE DETECTION SENSITIVITY AND OPTIMIZATION 



The vacuum radiation input noise b in is delta correlated and the thermal fluctuating force is approximated as white 
noise thus is also delta correlated, so we have the following correlation relations in frequency domain: 



(X in (uj)X in (uj')) = (Y in (uj)Y in (uj')) = -8(u) + bj'), 

(X in (u)Y in (u')) = -{Y in {u>)X in (oj')) = l*(w+w'), 
{F th (u)F th (u>')) = 2mjk B TS(u + Lj r ), 



(26a) 

(26b) 
(26c) 



from which we have 



7](ui) = 2m^k B T + - 



2mr(k B T 



2hy^G (k - iAf 



sm u — % cos ( 



2m"fkBT 



k — iuj + iA 
2H 2 kG 2 



2^/kG A sin 6 + (k — iu) cos 
{n-iAf +lo 2 \ 2 m 2 |< 2 -ui 2 -iiui\ 2 



| K -iw + iA| 2 8kG 2 
sin 9 + i cos 9 



H — (n + iA + iuj)- 

2 A sin 9 + {k + ilo) cos 



| A sin 9 + (k — iu) cos 9^ 
m(u}' m 2 — uj 2 + ij'uj) + c.c. 



(27) 



The first term is the thermal white noise. The second term proportional to G 2 comes the random back action force. 
The third term proportional to 1 /G 2 is the phase noise related to position measurement imprecision. 
In the vicinity of 9o = — arctan(«;/A), 9 = 9q + 69, and the optimal pump is approximated as 



which leads to an effective mechanical frequency 



I 2 _ 2 



2hAuj 2 a 2 /m,L 2 

t -2 T A2 



2k, 



(28) 



(29) 



The susceptibilities at DC reduce to 



XF 



\Xx - IXY\ 



1 A 

2hmuj 2 k 



(1 - 



(^+i)69 



(30a) 



(30b) 



In the limit 89 — > 0, the force susceptibility remains finite, while the optical susceptibilities go to zero. So the total 
sensitivity is given by 



I 4^ + A~ ) 



(31) 



